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We study the rotational properties of an attractively interacting Bose gas in a quadratic + quartic 
potential. The low-lying modes of both rotational ground state configurations, namely the vortex 
and the center of mass rotating states, are solved. The vortex excitation spectrum is positive for 
weak interactions but the lowest modes decrease rapidly to negative values when the interactions 
become stronger. The broken rotational symmetry involved in the center of mass rotating state 
induces the appearance of an extra zero-energy mode in the Bogoliubov spectrum. The excitations 
of the center of mass rotational state also demonstrate the coupling between the center of mass and 
relative motions. 



I. INTRODUCTION 

The appearance of quantized vortices offers a strong 
evidence of superfluidity in a rotating Bose-Einstein con- 
densed system. A classical example is the quantized cir- 
culation in superfluid ^He. More recently, the existence of 
these novel quantum states has also been demonstrated 
in dilute atomic gases 0, • It is found that above a cer- 
tain critical rotation frequency, atomic vapor with repul- 
sive interactions responds to rotation by forming a singly 
quantized vortex ^ As the rotation rate increases, 
more vortices appear and vortex lattices consisting hun- 
dreds of vortices may be detected 0, 0] . In addition to 
atomic boson gases, the experimental detection of vortex 
lattices in both sides of Feshbach resonance in quantum 
degenerate Fermi gas has just been reported 0. 

The response of a Bose condensed gas to an external 
rotation is expected to be different in the case of attrac- 
tive interactions. In a harmonic-oscillator potential, it 
is not the vortex state which is energetically favorable. 
Instead, the least-energy configuration for an attractive 
condensate with non-zero angular momentum is a state 
with a rotating center of mass (cm.) ^, 9,10]. However, 
this state is thermodynamically unstable, for the criti- 
cal rotation frequency to excite the cm. mode is equal 
to the trap frequency. The condensate set on rotation 
thus experiences the centrifugal force being larger than 
the confining force of the trap and eventually escapes the 
trap. 

Although a rotating attractive condensate is not stable 
in a harmonic potential, it can be stabilized by using a 
steeper confinement. A possible scheme is an anharmonic 
trap configuration in which a small quartic potential is 
superposed onto the harmonic trap. This potential is ex- 
perimentally realistic; it has already been implemented to 
study the fast rotation of repulsive *^Rb condensate ■ 
On the attractive side it is expected that the phase-space 
consists of states of broken cylindrical symmetry and, 
for very weak interactions, multiply quantized vortices 
[T^ Ha . Q, The states of broken cylindrical sym- 
metry involve a rotating center of mass and are closely 
connected with the cm. rotating state in a harmonic 



trap. Even though the correspondence is never exact, in 
this paper we use the term cm. rotating state for all the 
states with a broken cylindrical symmetry. 

The rich collection of rotational states of an attractive 
condensate in an anharmonic trap is worth further stud- 
ies. An especially interesting effect is the spontaneously 
broken cylindrical symmetry in the case of cm. rotating 
state. In addition, the investigation of dynamical and 
thermodynamical processes is a necessary extension to 
the work performed so far. This includes the study of 
the lack of superfluidity in attractive condensates in- 
dicating the fact that, if the external rotation is stopped, 
the condensate relaxes to the non-rotating ground state, 
independently of the initial rotational state. Recently, 
these topics have been studied in one dimension where 
even exact results have been obtained 17, 18, 19, 201 1. 

In the present work, we study the elementary excita- 
tions of rotating attractive Bose condensed gas confined 
in an anharmonic quartic -t- quadratic trap potential. We 
present our system and the relevant formalism in Sec. ITU 
The nature of the ground state excitation spectrum is 
discussed in Sec. IIIII and the results are summarized in 
SecHVl 



II. THE MODEL 

We assume a zero-temperature Bose gas with negative 
s-wave interaction. The anharmonic confining potential 
has the following form 
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which consists of the harmonic oscillator potential plus 
the radial quartic term described by the dimensionless 
parameter A. The oscillator length is defined as Oosc = 

Following the argumentation of Ref. jl^ . by ignoring 
the quantum fluctuations of the possible cm. motion, 
the ground state behavior of a rotated attractive Bose 
condensate can be approximated by the ordinary single- 
component Gross-Pitaevskii (GP) equation in the frame 
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of reference rotating with frequency ft. The GP equation 
for a stationary state is 



Ho-n-L + Uo\-9 (r) |2 ^^J (r) = /i* (r 



(2) 



where Hn — — T^—yS ^ y /y) is the Hamihonian for a sin- 
gle atom in the nonrotating trap. The interaction con- 
stant Uq is proportional to the scattering length by the 
relation Uq = Airh'^a/m. L = —ifi (r x V) is the angular 
momentum operator and /i is the chemical potential. The 
condensate wave function is normalized to the number 
of atoms N which is fixed here to = 1000. We also 
fix the mass m to the mass of atomic ^Li. For the an- 
gular frequencies of the harmonic oscillation we set the 
values cj = 27r X 30 Hz and = 27r x 180 Hz. Here, 
Jl is always pointing to the z-direction, which, together 
with the strong confinement in this direction, makes the 
problem effectively two-dimensional. 

The phase-diagram of the ground state solutions of 
Eq. is studied in Here, we extend our earlier 

work and study the spectrum of normal modes of the 
solutions of Eq. |(2Jl. These are obtained by solving the 
coupled Bogoliubov-de Gennes equations for the quasi- 
particles Uj and Wj, 



(3) 



where K ^ Hq - U ■ L + 2f/o|* (r) |^ " A* and Ej is the 
energy of the jth eigenmode. 

The calculation of the spectrum of the eigenmodes in 
our numerical procedure consists of first solving the low- 
est energy configuration of Eq. |(2Jl in a two-dimensional 
Cartesian grid by propagating the wave function in imag- 
inary time. The second step is to transfer the obtained 
wave function and chemical potential to Eqs. lO and 
numerically diagonalize these in order to solve Uj , Vj and 
E,. 



III. GROUND STATE EXCITATIONS 





FIG. 1: The lowest quasiparticle eigenenergies for ground 
state configurations calculated for different values of the scat- 
tering length a. In the upper graph (A — 0.15, il/tj = 1.35), 
and (A — 0.05, Q/uj = 1.15) in the lower graph. The dotted 
vertical line represents the interpolated value of the critical 
scattering length separating the m — 2 vortex ground state 
(right) and the cm. ground state (left). 



We analyze the numerical solutions of Eq. |j3Jl by fixing 
the pair of parameters (A, fi/w) and varying the inter- 
action strength Uq. We first choose (A = 0.15, ri/w = 
1.35). For these parameters the phase-space was ana- 
lyzed in Ref. In addition, we choose another pair 
(A = 0.05, ^/lo — 1.15) to be studied, for these val- 
ues are closer to the pure harmonic oscillator. The low- 
energy spectra are shown in Fig. ^ where we have varied 
the scattering length a to cover the phase-space portions 
of vortex and cm. states. In both cases, the ground 
state on the vortex side is doubly quantized. A vortex 
is only stable for weak interactions, and the increase in 
the magnitude of interaction strength leads to a phase 
transition at certain critical a = Uc- For these relatively 
strong interactions, the ground state is the one with bro- 
ken rotational symmetry. For very strong interactions 
(or high N) the system becomes unstable. 



A. Vortex ground state 

On the vortex side, the ground state as well as the ex- 
citations are cylindrically symmetric. Consequently, the 
Bogohubov amplitudes can be written as Uq cx e'('+™)^ 
with Vq oc Q^ii-™)^ ^ Here m (not to be confused with 
the mass) is the angular momentum quantum number of 
the condensate, which, in the case of a doubly quantized 
vortex is TO = 2. The excitation winding number q is 
the angular momentum relative to the condensate. En- 
ergy spectra in Fig. ^ show that the lowest excitations 
are the modes with q = I and q = —1. Especially, the 
energy eigenvalue of q = 1 mode is very low. Increasing 
a I decreases the energy of the mode quickly to negative 
values. The crossover determines the critical value Oc 
below which the vortex state is not stable. At a Qc 
the mode energy is real within the numerical precision. 
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FIG. 2; Density plots of diflerent wave functions related to 
the cm. state (upper left graph). The relevant parameters 
are A = 0.15, Q,/uj = 1.35 and a = —10.0 a.u. The upper right 
graph is the Goldstone mode <^g as defined in Eq. ||1J. The 
part of orthogonal to the cm. rotational wave function is 
denoted as x a^nd plotted in the lower left graph. In the lower 
right graph we plot the zero-energy Bogoliubov solution uq. 
The insets show the phase profiles of the wave-functions. 



However, in the unstable region the imaginary part of the 
negative mode rapidly increases as a function of the inter- 
action strength. This indicates the dynamical instability 
of the vortex state. We have checked this by performing 
time evolution GP simulations. Indeed, we have seen the 
breakup of the vortex wave function when \a\ is strong 
enough whereas for \a\ < Oc the vortex state is stable in 
time. 



B. Cm. ground state 

On the left hand side of the Fig. ^the excitations are 
calculated for the ground state solution which is the cm. 
rotational state. The ground state stability is ensured 
by the fact that the excitation energies are real and non- 
negative. We have also propagated the GP ground state 
solutions in real-time and we have not seen any sign of a 
dynamical instability. 

As seen in Fig. ^ there is an extra zero-energy mode on 
the cm. side of the phase-space (denoted by G). Due to 
Goldstone's theorem j^, the spontaneously broken sym- 
metry implies the existence of a zero energy mode. Gen- 
erally, the mean-field approximation breaks the global 
U(l) symmetry of the system [S^. However, in the case 
of cm. rotational state, the cylindrical 0(2) symme- 
try is also spontaneously broken. For a mean-field wave 
function (pc with a particular broken symmetry ipQ, the 
Goldstone mode is ^'''"^^'^"'^ In this case, the Gold- 
stone mode corresponding to the rotational asymmetry 



*e (r) = 



5^- (r) 
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In Fig. 121 we show an example of a cm. rotational wave 
function and the corresponding Goldstone mode (r). 
The density profile of $e (r) is slightly more elongated, 
but otherwise it is rather similar to I^P (r)|^. The struc- 
ture of ^0 (r) may be confusing at first sight, but there 
is an extra twist in the problem. Namely, if more than 
one symmetry is spontaneously broken, the correspond- 
ing Goldstone modes are not necessarily orthogonal. This 
means that the mode in Eq. is generally a superposi- 
tion of the cm. rotational wave function ^E* and a wave 
function x orthogonal to 5*. In this context, we write 



$e (r) = /3* (r) + x (r) 



(5) 



where /3 is a complex constant. Because of orthogonality, 
both X ^-iid [3 can be easily solved. In the example of Fig. 
12 we present the function x i^i the left lower graph. The 
extra zero-energy Bogoliubov solution for the same wave 
function can also be written as a superposition of the 
two orthogonal modes. Also, the mode amplitude uq 
has a component orthogonal to the cm. rotational wave 
function. This is plotted in the right lower graph in Fig. 
121 The identical structure (within numerical precision) 
with the fimction x is clearly visible. 

The differential operation performed on the conden- 
sate wave function in Eq. 10} is just iL. Although the 
function x is generally small compared to j3'^ it is clear 
that 5" is not an angular momentum eigenstate. On the 
other hand, the true many-body cm. rotational state in 
a harmonic trap @ is an eigenstate of L. The contradic- 
tion is a side effect of the broken symmetry, which, we 
believe, in turn is a consequence of the used GP approx- 
imation. This is supported by the quasi-one-dimensional 
analog 0, 0, |23| , where the symmetry is broken inde- 
pendent of the details of the trap potential which forms 
the torus. 

For actual, positive energy excitations the general be- 
havior is that the energy eigenvalues increase with the 
interaction strength. However, interpretations of individ- 
ual modes should be done with care. In one-dimensional 
torus, the lowest mode is a breathing mode [23|. Gener- 
ally, this is not the case in 2D because the anharmonic 
term in the trap potential couples the cm. and rela- 
tive motions. A change in the mean size of the cloud 
affects the effective potential seen by the cm. Con- 
sequently, a time dependent wave function correspond- 
ing to an excited mode is expected to show some cm. 
motion also. For the lowest excitation (denoted by A 
in Fig. the coupling can be seen in Fig. O where 
we have assumed an excited wave function of the form 



^ (r,i) = 5- (r) -f UAe 
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As a rough 



description of the oscillation, at t/T = the elongated 
wave function is on azimuthal motion. The cm. mo- 
tion then begins to decelerate and excitational breathing 
begins to compress the cloud. The compression rate is 
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FIG. 3: A cm. rotational ground state excited with the low- 
est true excitation (A) obtained from Bogoliubov-de Gennes 
equations. Here A = 0.15 and a = —5.0 a.u. and T ~ 78 l/uj 
is the period of oscillation. 



highest at the classical turning point of the cm. On the 
backward cm. motion the cloud retains the squeezed 
shape (see t/T = 0.35 plot in Fig. When approach- 
ing the other turning point the cm. motion slows down 
and the cloud begins to enlarge again. The coupling be- 
tween cm. and shape oscillations is not characteristic 
only to this mode but also to all other modes we have 
investigated. 

The investigation of the mode profiles for Uj and vj 
also reveals that the broken rotational symmetry is not 
as visible for the high energy modes as it is for the low- 
lying ones. In the computations we have performed, the 
mode B corresponding to the q — 4: mode on the vortex 
side carries only a slight asymmetry on the cm. side 



whereas the lowest excitation is highly localized. 



IV. CONCLUSIONS 

In summary, we have studied the Bogoliubov-de 
Gennes excitation spectrum of a rotating attractively in- 
teracting Bose-Einstein condensate trapped in an anhar- 
monic potential. The stability of the (multiply) quan- 
tized vortex state is ensured by the positive energy of 
the lowest excitations. However, the eigenmodes q = 1 
and q = —1 become unstable at relatively weak values of 
the scattering length. These modes being unstable the 
true ground state becomes the state with the rotating 
center of mass. For the cm. state the rotational sym- 
metry is broken and the Bogoliubov spectrum reveals an 
extra zero-energy eigenmode. This Goldstone mode is 
the rotation mode of the cm. state. Finally, by exciting 
the cm. state with lowest positive energy mode, we have 
illustrated the coupling between the cm. and relative 
motions in anharmonic potential. 
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